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ABSTRACT 

Qh| In a recent article, Faraoni proposed an alternative procedure to solve the 

O ' Friedman-Lemaitre-Robertson- Walker (FLRW) cosmological equations. The 

basic result of that paper was obtained long ago through a different approach, 
! which seems to be little known and deserves closer attention due to its peda- 

gogical interest. The broad importance of this method is readily recognized 
^ by examining some additional cases not considered by Faraoni. Its instructive 

■ potential for introductory courses on cosmology has been positively verified 



by the author in many lectures during the last decade. 



Recently, Faraoni|]I| rediscussed the cosmological solutions for a relativis- 
tic simple fluid in the framework of FLRW models. In this note we present a 
variant of the method proposed by Faraoni, which seems to be more general 
and quite interesting from a pedagogical viewpoint. Although little known, 
this approach was proposed long ago [0], and its pedagogical efficiency has 
also been verified many times in the last decade. In what follows, after a 
comparison with the results presented by Faraoni, we stress the generality of 
this procedure by discussing a class of solutions with a cosmological constant. 
For the sake of completness and further reference, as well as to consider some 
subtleties not addressed in Faraoni's work, we first summarize the basics of 
this problem. 

The spacetime metric takes the following form (c = 1) 

ds' = dt" - a{tf { + r^de^ + r^stn^edA (1) 
— Ar^ J 

where a{t) is the scale factor, and K = 0, ±1 is the curvature parameter. In 

such a background, the Einstein field equations for a relativistic simple fiuid 

can be written as [1-5] 

^ d^ K 
MGp = 3— + 3— , 2 

d K 

8nGp = -2 ^ - ^ , 3 

where p and p are the energy density and pressure, respectively. In this 
system there are three unknown quantities, namely a{t), p{t) and p{t) and 
only two independent equations. Thus, in order to solve it, an additional 
constraint is required. In the cosmological framework, it is usually assumed 
that the matter content obeys the barotropic 7-law equation of state 
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p=(7-l)p, 7e[0,2] . (4) 

This happens because such an expression provides a rather simple description 
for several cases of physical interest 0: (i) dust or pressureless matter (7 = 
1), (ii) radiation (7 = 4/3), (iii) vacuum (7 = 0), and (iv) Zeldovich's stiff 
matter (7 = 2). Naturally, the above assumption does not imply that the 
7 — parameter remains constant in the course of the cosmological evolution. 
In the standard hot big bang description, the Universe underwent a transition 
from a primordial radiation phase (7 = 4/3) to a matter-dominated epoch 
(7 = 1) during the late stages of its evolution. 

The cosmic dynamics is determined by combining the above set of equa- 
tions. In principle, the corresponding dynamic behavior must be heavily 
dependent on the choice of the two free parameters: (i) the curvature pa- 
rameter K, and (ii) the equation of state parameter 7. As one may check, 
the evolution of the scale function is driven by the second order differential 
equation |1|, |^ 

ad + Ad^ + AK = , (5) 
where the parameter A (c in Faraoni's notation) is a function of 7 

. (a) 

The first integral of (|) is 

^'=(-T-K , (7) 



a 



where the convenient integration constant, Cq^^, requires that all curves de- 
scribed by the general solution share the point a = at some instant of the 



2 



cosmological time. Furthermore, the derivative a(ao) = VI - fsT is the same 
for all values of 7, and indicates the existence of an extremum for closed 
models. We see from (^) and @ that this extremum is a maximum if 7 > |, 
a minimum if 7 < |, or a stationary point if 7 = | 0. In addition, equation 
(ID tells us that the energy density 



3 /a„\37 



(8) 



is positive definite, as should be expected for a physical fluid. Note also that 
for K = the first integral (|^ has the following solution 

a(t) = a„[l + (A + l)(t-to)/ao]^, (9) 

or equivalently, 

a{t) = ao l + -^( -) , (10) 

where to is an arbitrary time scale. In particular, by choosing to = 2ao/3'y one 
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recovers the restricted form of the flat solutions, namely, a(t) = ao{t/to)^^ ■ 
In addition, in the limit 7 — > 0, the above expression reduces to a{t) ~ e^°*, 
which is the de Sitter fiat solution {Ho = a~^). 

At this point we remark that the aim of Faraoni's paper is to obtain the 
general solution of the first integral (^, or equivalently of @. This can be 
accomplished by using the conformal time 1], instead of the cosmological or 
physical time {dt = a{ri)dri). In this case, the metric (|1]) and equation (|^) 
takes the forms below 

ds^ = a{r]f (^df]^ - Y^J^ - - rhin'9d4>^^ (11) 

aa" + (A - l)a'^ + AKa^ = , (12) 
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where the prime denotes derivative with respect to conformal time. Now, 
instead of using Faraoni's transformation, u = a' /a, which leads to Ricatti's 
differential equation (see section 3 of [|l|), we employ the auxiliary scale factor 



Zir]) = a^, tf A^O (13) 
Z{r]) = Ina, if A = , (14) 



to obtain, respectively. 



Z" + KA'^Z = 0, if A^O (15) 

Z" = 0, if A = . (16) 

Equation (|l^) is identical to equation (3.7) in Faraoni's paper. Note that 
it reduces to (^) in the limiting case A = 0. Therefore, it can be considered 
the general equation of motion for any value of A. As remarked before, it 
depends strongly on the pair of parameters (7, if). The physical meaning 
of (|T5|) is apparent: It describes the classical motion of a particle subject to 
a linear force. This force is of restoring or repulsive type depending on the 
sign of the curvature parameter. 

Closed models (K=l) are, for any value of A 7^ 0, analogous to harmonic 
simple oscillators (HSO). Therefore, the cosmic dynamics in this case is sim- 
ilar to a spring-mass system where the spring constant is determined by the 
7-parameter defining the equation of state obeyed by the cosmic fluid (see 
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Eq-(@))- It is worth noticing that for positive values of A, that is, 7 > 2/3, 
this general and intermittent oscillatory motion between both singularities 
("big-bang" and "big-crunch") reinforces the connection with the idea of a 
pulsating Universe. From a pedagogical viewpoint, these harmonic solutions 
for closed models are more enlightening than the cycloidal parametric so- 
lutions usually presented in textbooks for particular values of 7 [3-5]. It is 
actually surprising that the dynamics of closed universes might be reduced to 
that of a simple harmonic oscillator, a key system for analytical and algebraic 
calculations in many disparate branches of physics. It also leads naturally to 
the question whether our Universe is the largest clock, that is, a global os- 
cillator. Indeed, although the current literature (including many textbooks) 
admits only a unique incomplete cycle to describe the observed Universe, the 
existence of pulsating solutions fascinate many philosophers and cosmologists 
since they are associated with the old concept of an ethernal return. 

It is also worth mentioning that if the Universe is spatially flat [K = 0), 
equation (^) implies that the system behaves like a free particle, and the 
same happens if A = 0. In the later case, this effective free particle behavior 
holds regardless of the curvature parameter. For hyperbolic spacetimes {K = 
— 1), the system behaves like an "anti-oscillator", that is, a particle subject 
to a repulsive force proportional to the distance. 

Now, recalling the elementary mathematical identities: 



lim 



smax 



(17) 
(18) 



= X 



a 




the unified solution of (|T5D can be written as 




(19) 
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where Zo and 5 are integration constants. Note also that using the transfor- 
mation (|T3p, the first integral ([7p may be translated into the energy equation 

2 2 2 ^ ^ 

which determines Zq = ao^- 

By choosing 5 = 0, the general parametric solution relating the scale 
factor and the cosmological time, dt = a{ri)dri, is given by 



a{r])=ao[ TtT^] ' (21) 



K 



, , , , sinV K\A\ri\'^ . 
t[ri) = tto [ 1= + constant . (22) 



K 

From identity (18) these solutions for elliptic and hyperbolic models read: 
K = 1 



a{r]) = ao{sin\A\r])'^ (23) 
t{ri) = ttoj {sin\A\ri)'^ + constant (24) 



K = -1 



a{r]) = Go {sinh\A\r]) ^ (25) 
t{r]) = ttoj {sinh\A\ri)^ + constant (26) 

where the integrals (|2^ ) and (|26| ) for t{ri) may be represented in terms of 
hypergeometric Gaussian functions |@|. As one may check using (17), the 
integral for i^' = is trivial, and the parametric solutions may readily be 
inverted to give the scale factor as a function of the cosmological time. 

It is also interesting to show that the above method based on the trans- 
forming equations (p!3| ) and (|T^) is also convenient when new ingredients are 
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considered such as the presence of a cosmological A-term In this case, 
the EFE equations can be written as 





K 






+ 3- 




d 


a" 


K 


2- - 
a 


' a2 ~ 





(27) 
(28) 

where A is the cosmological constant. Now, using (|), one may see that the 
scale factor satisfies the generalized FLRW equation 

ad + Ad^ + AK~ \a{A + l)Aa^ = . (29) 

Following the same steps of the case with no A, instead of (|T2D, we obtain 
in the conformal time 

aa" + (A - l)a'^ + AKa^ - -A{A + l)Aa^ = , (30) 

3 

and transforming to the auxiliary scale factor through (|T3|) results 

Z" + KA^Z - -A{A + l)AZ^ = . (31) 
3 

As expected, if A = 0, equations (|2|), (jSg) and reduce to (|), (|12]) and 
(|15D, respectively. 

Physically, the equation of motion ( ^T]) means that closed universes with 
cosmological constant evolve like anharmonic or non-linear oscillators [p!0| 



The anharmonic contribution to the HSO is proportional to the cosmological 
A-term, and its power index depends uniquely on the equation of state 7- 
parameter. Note that the method provides an exact non-linear description. 
In other words, the anharmonic term does not appear due to the retention 
of higher order terms in some particular perturbative scheme. 



7 



Finally, we remark that the unidimensional equation of motion (^), as well 
as its generalized version with a cosmological constant are endowed with an 
obvious mechanical analogy which can be described within the classical La- 
grangian formalism. It has been shown that a particle subject to a potential 
V{a) = aa~"' — |Aa^, where a and n are constants, may satisfy the equation 
of motion (pQ]). More precisely, this happens if n = 2 A, and the energy 
is proportional to the curvature, 2E = —rriK, where m is the mass of the 
test particle 0. Naturally, the same method of solution, shown here for the 
relativistic case, may also be applied to this classical Lagrangian approach. 
The author has explored both approaches in many introductory courses on 
cosmology, and verified the nice pedagogical impact on the students. 
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